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Short strings and gluon propagator in the infrared region. 
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We discuss how infrared region influence on short distance physics via new object, called "short string". This 
object exists in confining theories and violates the operator product expansion. Most analytical results are 
obtained for the dual Abelian Higgs theory, while phenomenological arguments are given for QCD. 



1. Introduction 

Below we describe effects produced by the re- 
cently introduced object, short string. Short 
strings ||l| are responsible for non-standard power 
corrections in QCD, which are important in the 
ultraviolet region. On the other hand short 
strings exist due to confinement and in this sense 
they are induced by large distance physics. Short 
strings are not point-like objects and violate such 
fundamental principle as operator product expan- 
sion (OPE). 

It is interesting to find the short string effects 
in QCD, but confinement in QCD can not be de- 
scribed from first principles and at the beginning 
(in Sects. 2-4) we consider the dual Abelian Higgs 
model (AHM). In this model electric charges are 
confined at the classical level, and it is simple 
enough to study both large and small distances. 
Moreover, due to the Abelian dominance [||, the 
dual AHM can be considered as the effective in- 
frared theory for QCD. At small distances dual 
AHM is not related to QCD, and we consider it 
to study general properties of short strings. We 
show that short strings give rise to special singu- 
larities in the photon propagator and produce the 
linear correction, agr, to the Coulomb potential, 
a/r, at small distances. 

In Sect. 5 we discuss numerical (lattice calcula- 
tions) and phenomenological facts supporting the 



existence of short strings in QCD. 

2. Short strings in classical equations of 
motion for AHM 

Consider the following problem: what are cor- 
rections to Coulomb potential at small distances 
in QCD? If we neglect the trivial perturbative de- 
pendence of the bare coupling on r, we have in the 
standard approach (perturbation theory -I- OPE): 



Vir) 



-a/r + c ■ r'^ , and 7 = 2 or 3 depending 



on the parameters of the system (see |3| 



)• 



Consider now the dual AHM with the action: 



(1) 



where B is the dual gauge field, $ is the Higgs 
field carrying the magnetic charge g. The Higgs 
field and the gauge field are massive: mj^ = 
2\rf,my = g^rf. The classical potential be- 
tween +e and — e charges (analogue of the quark 
- antiquark pair) at large distances is: 
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*Material based on talks by F.V. Gubarev and M.I. Po- 
likarpov at Lattice '99, Pisa, Italy. 



The confining behavior is due to the dual 
Abrikosov - Nielsen -Olesen (ANO) string which 
is formed between -\-e and — e charges. At small 
distances the potential is Coulomb like: V{r) = 
—a/r, r — * 0. We are interested in the classi- 
cal corrections to the Coulomb potential at r ^ 
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y. One has to impose the boundary 
condition, that is the vanishing of the scalar field 
along a mathematically thin line connecting the 
magnetic charges. This topological condition was 
formulated in Ref. |^ and implied in numerical 
studies of the problem (see, e.g., Q). 

This infinitely thin line is nothing else but the 
Dirac string connecting the external charges. By 
definition the energy of the Dirac string is zero. 
It is easy to realize that the Dirac string cannot 
coexist with $ 7^ 0. Indeed, if the Dirac string 
would be embedded into a vacuum with ($) ^ 
then its energy would again jump to infinity since 
there is the term 1/2 in the action and 

— ^ 00 for a Dirac string. Therefore, even in 
the limit r — > there is a deep well in the profile 
of the Higgs field $. This might cost energy which 
is linear with r even at small r. 

Now we are coming to the next question, 
whether this mathematically thin line realizes as 
a short physical string. By the physical string we 
understand a stringy piece in the potential, CTo ■ r 
at small r. In other words, we are going to see 
whether the stringy boundary condition implies 
a stringy potential. To get the answer one solves 
the classical equations of motion, the result is m : 



lim V{r) 



r 



+ (7o-r + 



(3) 



with a non- vanishing slope a^. 

The slope ctq depends smoothly on the value of 
my/rriH- For example for my — ttih the slope of 
the potential at r — > is the same as at r ^ cx). 
That is, within error bars: 



(4) 



where a^o determines the value of the potential 
at large (^. 

3. Gauge boson propagator in the confin- 
ing theory 

The simplest "phenomenological" theory of the 
quark confinement is formulated in terms of the 
gluon propagator. Namely, one argues that 
if the gluon propagator has a double pole in the 
infrared region then the quark-antiquark poten- 
tial gets a linear attractive contribution at large 



distances and this implies the quark confinement. 
On the other hand, Gribov and Zwanzigcr ar- 
gued that elimination of the gauge copies leads to 
a counter-intuitive vanishing of the gluon prop- 
agator in the infrared region. Some analytical 
studies [|l^ confirm this behavior, while the oth- 
ers [ pl predict a singular gluon propagator in 
the infrared limit. Recent numerical simulations 
in lattice gluodynamics ||l2| show that the gluon 
propagator seems to be finite at small momenta. 

Below we consider the gauge boson propaga- 
tor in dual AHM, that is the abelian analogue 
of the gluon propagator. To simplify calculations 
we consider the London limit, that is A — > cx) 
in (|l|), and mu — *■ 00. The final result remains 
unchanged if we relax the condition on the cou- 
pling A. To evaluate the propagator we utilize 
the Zwanziger local field theory of electrically 
and magnetically charged particles . The the- 
ory contains two vectors potentials, namely the 
gauge field A^j^{x) and the dual gauge field Bf^{x), 
which interact covariantly with electric and mag- 
netic currents, respectively. The corresponding 
Lagrangian is: 



(5) 



.zw {A,B)+ ie il [x) Af, {x) + 

+igjS{x)B^{x), 

where e (g) stands for the electric (magnetic) 
charge, and (j^) is the electric (magnetic) ex- 
ternal current. The Zwanziger Lagrangian Czw 
is rai: 
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{n-[dAA]) 
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in-[dAB]f 



+ ^{n ■ [d A A]){n ■ [d A Bf) 



--in-[dAB])in.[dAAf), 



where \A A B 



fJ.V 



n^{A A B)^^ ^ (G)^^ = ie^^Ap 



Ap 



(6) 



[A A 
. In- 



spite the Lagrangian (j|) contains two gauge fields, 
A^ and B^, there is only one physical particle (a 

The 



massless photon) in the spectrum [ 13[ 
Lagrangian ^ depends on an arbitrary constant 
unit vector Ufj,, nj^ = 1, while physical observables 
are insensitive to the direction of n provided the 
Dirac quantization condition, e-g — 27rm, m G Z , 
is satisfied [Ol. 
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The theory of condensed monopoles interacting 
with an external electric source j"^ (quark cur- 
rent) is described by the following partition func- 
tion 

Z[f] ^ J VAVBV^ exp|-y" d'^x(^Czu:iA, B) 

+ + igB)<S>\' + ~ r^^f ~ lej^A) } . 

We work in the Euclidean space and consider the 
London limit, A ^ oo. The gauge field propaga- 
tor is defined as: 
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5j'L{x)5jl(y) Z[0] 



(7) 



To specify the propagator completely we choose 
the axial gauge n^Af^ — 0. The standard param- 
eterization of the propagator in this gauge is js) : 



D 



axial 



+ 



{p-n) 



{P -n) (jp. n) 

{6^..-n^n,)G{p'). (8) 



Using methods of Refs. (l^ the functions F(p^) 
and G{p^) can be expressed through the ANO 
string-string correlator : 



G{p') = - 
where 
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m\,D^{p'^ 



p2 D^ip') 



(p2 + m.y){p ■ n) 



(9) 
(10) 

(11) 



(fc2 + my){k ■ n) 



and the correlation function of ANO strings 
< Tifj,a{p)'^vi3{—k) > is formally defined in [^6). 

The propagator (^ has singularities not only 
in the p'^ plane but in the variable {p ■ n) as well. 
The singular in (ja-n) terms appear not only in the 
longitudinal structures but also in front of {S^i, — 
n^rii,). It turns out that the singular structure is 
due to the Dirac strings, which are directed along 



vector entering the Zwanziger Lagrangian (^ . 
Note that Dirac strings become short strings or 
ANO strings in the dual AHM. 

Although expression (p[]Tl|) for the gluon prop- 
agator is exact, it contains an unknown function 
-D^(p^). The infrared behavior of this function 
however can be qualitatively understood on phys- 
ical grounds. Indeed, consider the string part 
( |Tl| ) of the gluon propagator (^). This part corre- 
sponds to the following process: the quark emits 
a gauge boson A which transforms to a dual gauge 
boson B via coupling < A ■ B > which exists in 
the Zwanziger Lagrangian (|^). The dual gauge 
boson transforms back to the gauge boson after 
scattering on a closed ANO string world sheet 
E. The intermediate string state is described by 
the function D^(p^) and this state can be consid- 
ered as a glueball state with the photon quantum 
numbers 1~. 

The behavior of the function D^{p^) in the 
infrared region, p 
D^ip')^^^ + 
less parameter and Mg\ is the mass of the lowest 
1^ glueball. The dots denote the contributions of 
heavier states. Thus, according to eqs. (^,^0|) in 
the infrared limit, p 



0, can be estimated as 
where C is a dimension- 
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C 



F''\p') = — + T7T 
Gi^(p2) = -1 



(12) 



Note that neither double, l/(p^)^, nor ordinary, 
1/p^, pole are exhibited by the gluon propagator 



- |12[ ) in the infrared limit. There are gauge de- 
pendent l/(p • n) singularities which reflect the 
presence of Dirac strings. Therefore our gauge 
boson propagator does not vanish in the infrared 
either. In view of this, the vanishing of the propa- 
gator at p ^ predicted in |^ seems to be specific 
for the special gauge choice in QCD which leads 
to a nontrivial Fundamental Modular Region. 

4. Relation of potential and propagator, 
violation of the operator product expan- 
sion. 

As it is mentioned above, the existence of short 
strings is manifested also through breaking of the 
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standard OPE. The physical reason of this break- 
ing is that strings have extra dimensionahty com- 
pared to particles, they are not point-Hke objects. 
The another manifestation of this problem is the 
violation of the standard relation between the po- 
tential and the propagator. Let us try to recon- 
struct the classical Coulomb potential between 
(heavy) electric charge-anticharge and between 
magnetic charge-anticharge in the dual AHM for 
r <C rriy^, ra^ . For magnetic charges the po- 
tential can be reconstructed in the standard way 
from the propagator Dfii^{q^) in the momentum 
space, 

V{r)^ J d\ e^-' b^^ici^) . (13) 

In order to have the standard notations we use 
here q in the Minkovsky region, = —p^ (cf. eq. 
(^ - (pl|)). For q in Euclidean region and much 
larger than the mass parameters, the propaga- 
tor D^^{q^) can be evaluated by using the OPE. 
Restriction to the classical approximation implies 
that loop contributions are not included. How- 
ever, vacuum fields which are soft on the scale of 
can be consistently accounted for in this way 
(for a review of see, e.g., iQ). Modulus longitu- 
dinal terms, we have: 



2 

If the potential energy is given by the Fourier 
transform of ( p^ then its dependence on is ex- 
plicit. Thus it is impossible to get the Coulomb 
potential by the Fourier transform (^3|), we ad- 
dressed this problem in Ref . . 

Note that Eq. (|l^) immediately implies that 
the standard OPE does not work any longer on 
the level of q~^ corrections. Indeed, choosing q^ 
large and negative does not guarantee now that 
the my correction is small since the factor [qnf 
in the denominator may become zero. 

The reason for the breaking of the standard 
OPE is that even at short distances the dynamics 
of the short strings should be accounted for ex- 
plicitly. In particular, in the classical approxima- 
tion the string lies along the straight line connect- 
ing the magnetic charges and affects the solution 
through the corresponding boundary condition, 
see above. The OPE is valid as far as the probe 
particles do not change the vacuum fields dras- 
tically and the unperturbed vacuum fields are a 
reasonable zero-order approximation. In our case, 
however, the Higgs field is brought down to zero 
along the string and this is a nonpcrturbative ef- 
fect. 



+ = (14) 

q'^ (7^ / q — riiy 

Thus, one uses first the general OPE assuming 
\q^\ » g^'^^ then substitutes the vacuum ex- 
pectation of the Higgs field $ and upon summa- 
tion of the whole series of the power corrections 
reproduces the propagator of a massive particle 
(< >= rf, my — g^if). This propagator ac- 
cording to (^^ corresponds to Yukawa potential, 
or to Coulomb potential at small distances. 

The photon exchange between electric charges 
(analogues of quarks) corresponds to propagator 
(H). At the classical level we should neglect the 
contribution of virtual strings, D^ijP') ~ 0, and: 



5. Short strings in QCD 

Now we briefly review the possible manifesta- 
tions of short strings in numerical simulations of 
lattice gluodynamics. 

(i) The lattice simulation do not show 
any change in the slope of the full QQ poten- 
tial as the distances are changed from the largest 
to the smallest ones where Coulombic part be- 
comes dominant. Thus eq. (0) is presumably sat- 
isfied. Moreover, it is known from phenomeno- 
logical analysis and from the calculations on the 
lattice pO|-p^ that the realistic QCD corresponds 
to the case k « 1 where k, = mu/mv- It is re- 
markable that, as is mentioned above, the AHM 
in the classical approximation also results in the 
relation (^ for k « 1. 

(ii) Addressing the question of linear correc- 
tions at small distances one tries to explicitly sub- 
struct the perturbative contribution from the full 
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QQ potential. It appears that the remaining part 
is well parameterized by linear function and the 
extrapolation of the lattice data to very small dis- 
tances gives large value for (Tq, ctq ~ Bctoo p3| . 

(Hi) There exist lattice measurements of fine 
splitting of the QQ levels as a function of the 
heavy quark mass. The Voloshin-Leutwyler |p3| 
picture results in a particular pattern of the mass 
dependence of this splitting. Moreover, these pre- 
dictions are very different from the predictions 
based on adding a linear part to the Coulomb 
potential (BuchmuUer-Tye potential |2^). The 
results of most advanced measurements of this 
type are in favor of the linear (^^) correction 
to the potential at short distances. 

{iv) There is an evidence that the nonperturba- 
tive fluctuations on the lattice responsible for the 
confinement can be identified with the so called P 
vortices . If one measures the quark potential 
due only to these vortices, the numerical results 
seem to indicate that the slope of the potential is 
the same at large and at small distances, see [p7| . 

[v) The dimensional analysis show that the lin- 
ear correction at short distances to the poten- 
tial V{r) produced by short string corresponds 
to l/Q^ correction in OPE. The standard lead- 
ing correction is only l/Q*, which corresponds to 
the operator < TrF'^^j >. Numerical data for the 
expectation value of the plaquette minus pertur- 
bation theory contribution show 1/Q^ behavior 
(28|. 

At present the systematical errors in numerical 
calculations are large, especially the subtraction 
of perturbation theory contribution in items (ii) 
and (v) is not a well established procedure, never- 
theless it seems that there are several independent 
facts showing the existence of the short strings in 
lattice gluodynamics. 

Now we discuss two effects which can be man- 
ifestation of the short strings in QCD. 

(i) Analytical studies of the Bethe-Salpeter 
equation and comparison of the results with the 
charmonium spectrum data favor a nonvanishing 
linear correction to the potential at short dis- 
tances 

[ii) The most interesting manifestation of short 
strings is corrections to the standard OPE 
for current-current correlation function n((5^) 



H,^,^. It is impossible to calculate the coeffi- 
cient of 1/Q^ corrections, and in ref. [|o| it was 
suggested to simulate this correction by tahyonic 
gluon mass. The Yukawa potential with an imag- 
inary mass has the linear attractive piece at small 
distances as it is induced by short strings. The 
use of the gluon propagator with the imaginary 
gluon mass (m^ = —0.5 Gev^) unexpectively well 
explains the different behavior of 11 (Q^) in differ- 
ent channels. It is very important to perform the 
accurate calculations of various correlators 11 (Q^) 
on the lattice to check the model. There are the- 
oretical schemes which imitate the tachyonic 
gluon mass. 
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